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Numerical integration of the differential equations describing charge transport in plane-parallel 
ionization chambers allows to obtain charge profiles in addition to current-voltage curves. It is found that 
charge densities may exhibit maxima within a certain range of saturation. The effect is interpreted in terms 
of field modulation of carrier velocities. 
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1. Introduction 

The theory of the ionization chamber has been discussed in numerous papers since the 
turn of the century. Approximate solutions were given by Mie [l] 1 and Boag and Wilson [2]. The 
literature has been surveyed by Greening [3J. 

With the advent of modern computational methods precise numerical solutions can be 
obtained of many of the differential equations which until now defied attempts to obtain such 
solutions. 

By these methods, Armstrong-Tate [4] have solved the problem of the charge transport in 
parallel-plate ionization chamber. These authors calculated current-voltage curves and field 
profiles and compared their computer solution with the approximation given by previous authors. 
In this paper we use a computer method to calculate charge profiles. Results indicate novel 
features which so far have not been commented on. Our method of computation is similar to 
that of Armstrong-Tate. The use of different variables leads, however, to more advantageous 
conditions. Therefore, we have included a summary description of the calculus. 



AMSiMOS ISubjecl Classification .651)30 

•Current address. Stern wartenlr. 55, 1180 Wien. Austria. 

Figure* in bracket* indicate the literature references at the end of this paper. 

65 

571-453 O - 75 - 5 



2. Basic Equations 

We limit ourselves to steady-state conditions. The terms containing time disappear and one 
obtains the following well-known system of equations (for definition of symbols see table 1). 



edE'/dx' = p+ — p'_ 

di'+ldx' = - (a'le)p' + p'- + q 

- di'Jdx = - (a'/e)p'_ p' + + q 

*+ = ^+ P+ E' 

i'_ = pi pi E' 

E' = - dV'Idx' 



(1) 



These equations are completed by the usual boundary conditions: 



( J d ) E\x')dx' = V (2) 

pi (0) = p'_ (d) = 



By the introduction of the reduced, nondimensional quantities, (also defined in table 1) eqs 
(1) and (2) are transformed into eqs (3) and (4) which contain only 2 numerical parameters. Our 
equations are different from those of reference [4] because we use different reduction factors. 



dE/dx = p + - p_ 

di+/dx = - ap + p_ + 1 

- di_/dx = - ap + p_ + 1 

i + = p+ E 

/_ = Kp_ E 

E = - dV/dx 



(3) 



and 



( J d ) Edx = V 

P+ (0) = p_ (1) = < 
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(4) 



The eqs (3) may be transformed into 

dpJdV -- [p, (p 4 - p_) + ap 4 p 1|//: 2 

dEldV = - (p 4 - p.)IE 

Introducing as a new variable u = p+ — p_, one obtains 



(5) 



P_ = [(//£) - «]/(] + k) 
P + = [w + UkE\kI{\ + k) 

+ i r ^p+ i«i 

K L Jx k/:' 2 J 



(6) 



du k + l r Jp + 

Sc ~" ^r - L sr k£ 2 j (7) 



These relations can be simplified furthermore by introducing the variables E* — E/i and K* 
V/i 2 . Substitution of the eqs (6) and (7) into eq (5) allows to eliminate i and gives 



du 
dV* 


K+l 

kE* 2 


[ (1+K) 2 (Q ' K l) 


ua (k— 1) 
(k+1) 2 £> 






dE*/dV* = 


= - u/E* 



(k+1) 2 £> 2 



(8) 



3. Numerical Integration 

The variable // must have a zero at a point < x < 1; thus u(x ) = 0. Since p_(l) = 0, U 
increases monotonously from u = at x = x to « = p + at a: = 1. Thus i = p + (l)E(l) = M(i)£(l). 
With £* = £/i this gives u(l)E*{l) = 1. 

Equations (8) are now integrated for increasing values of u from u = to a maximum value 
which satisfies the relation wis* = 1 and corresponds to X = 1. Between a = .r and X = 0, u 
decreases from to a minimum value u(0) - - p_(0). The same argument used above gives 
Kp_(0)E(0) = i or - ku(0)E*(0) = 1. The integration terminates when the product uE* satisfies the 
relation - uE* = \/k which corresponds to x = 0. 

The value of the potential V is defined by postulating V*(x ) = 0. 

The integration is performed in terms of the positive parameter E*(x ) = E%. For E —* 
(linear region) one has p + = p_ and di/dx = 0. This gives p+p_ = 1/a and / = (p + + Kp_)E = (1 + 
K)Eja m or E* = oj 1/2 /(i + k). This is the lower limit of £>. 

This initial condition in terms of a known quantity assures that the system (8) always has a 
solution. Armstrong-Tate [4] use an initial condition which contains an assumed value of the 
field strength at one electrode. The equations of the problem admit however a solution only for 
values of the field higher than a certain minimum value. Since this is not known in advance, one 
is not sure to obtain a solution. 
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Therefore we have the boundary conditions 

u(x ) = 

< x =£= x : - (1/k) ^ uE* < 

x < x ^ 1 : < uE* ^ 1 

E*(*q) = E*m* « 1/2 /(l + k) < E* m 

V*(x ) = 

The current is obtained from the relation dV/dx = — E, giving 

dx = - i dV*/E* 
Integration over x for < x < 1 gives 

V*(l) 

i = - 1/ J dV*/E 



(9) 



(10) 



V*(0) 

(11) 



Since K* = V/i 2 , the applied voltage V = V(0) - V(\) is given by 

V = i 2 [V*(0) - K*(l)] 



(12) 



V*(0) and V*(l) are, respectively, the values of K* which are obtained when u has reached its 
extreme values given by the eqs (9). 



4. Current— Voltage Curves 

For air at NTP, Scott and Greening [5] have given /jl' + = 1.59 X 10~ 4 m 2 /T sec, /jl'_ = 2.21 X 
10" 4 m 2 /F sec, and a ratio ''recombination coefficient to elementary charge" = 4790 cm 3 /esu sec 
which in our notation give a reduced recombination coefficient a = 0.799 and k - 1.39. We have 
carried out the integration with these values. Results are represented in figure 1, curve A, which 
gives i\i % versus {i/i^/V 2 . 

In recent measurements Ritz and Attix [6] found a value /jl' = 1.89 X 10 4 m 2 IV sec for the 
mobilities of both positive and negative ions; their value of the ratio "recombination coefficient 
to elementary charge" was 4717 cm 3 /sec esu which in our notation gives a = 0.662. When the 
carrier mobilities are equal one has k = 1 and eqs (8) are considerably simplified, containing 
then only the single material parameter a. In addition the system is symmetrical with regard to 
the plane x = 1/2. Boundary conditions are also simplified. Figure 1, curve B, gives the 
computer solution for this case. In the saturation region curves A and B coincide. Even beyond 
this region, differences remain small. We have also carried out a series of integration for k - 1 
and values of a between 0.1 and 0.9. Results are given in figure 2. 
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Fi<;i KK 1. Saturation curves : A ) computer solution for p' + = 1 .59x10 4 m 2 / V sec; /j.' = 2.21 x JO 4 m 2 / \ sec ; a = 0.799. 
B) Computer solution for a = 0.662 and p' + = p.' = 1 .89x10 m \V sec. Circles and triangles indicate values actually 
calculated . 
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FiGt'RK 2. Saturation curves assuming equal mobilities ; computer solutions. Values of a: 0.1 (A), 0.3(B), 0.5 (C), 0.7 
(D),0.9 (E). C/i' + = /nl= 1 .S9xl()~ 4 m 2 l V sec) 



5. Charge Profiles 

Once i is known we can perform a second integration of eqs (8) and (10), obtaining the 
functions u(x), E*(x), and V* (x) . Using eqs (6) one also obtains the variation of positive and 
negative space charge densities with x. 



5.1. Equal Carrier Mobilities 

Figure 3 shows the positive space charge density as a function of x for equal mobilities and 
a = 0.662. The values of the negative charge density are symmetrical to the given curve with 
regard to x = 1/2. Different curves refer to different saturation ratios i/i s . For low and high 
values of this ratio the charge density is a monotonously increasing function of x but in the 
intermediate region a maximum is found to occur. Physically this can be understood as follows: 
as one moves from anode to cathode the density of positive charge is expected to increase due 
to the cumulative flow of positive charge carriers moving out from the anode region. The 
velocity of the carriers is, however, modulated by the field. In the center the field is weakened 
and the velocity reduced. In the cathode region the field is increased and the carriers sped up. 
Flow lines are "spread out." With increasing field strength the latter effect predominates and 
the positive charge density falls off. "Carriers are sucked into the cathode. " In the anode region 
the cumulative effect of carrier movement always predominates. Therefore, the curves have 
never more than one maximum. When approaching saturation the field becomes uniform and the 
effect disappears again. 

This view is confirmed if one solves the basic equations for a = 0, e.g., neglecting 
recombination; in this case one has always i/i s = 1. Boag and Wilson obtained approximate 
analytical solution, assuming that the mobility of the negative carriers is very great compared to 
that of the positive carriers. A precise analytical solution can, however, be obtained also in the 
general case where the mobilities have arbitrary values or are equal. We shall omit here the 
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Figure 3. Charge profiles positive space charge density versus depth . Computer solution for equal mobilities and a 
0.662. A )W S = 0.6,B;i/i s = 0.74,C)ili s =0.99 Q*; = (jl'_ = 1.89x10 4 m 2 / V sec) 



<? + 




FlGI KK 4. Charge profiles without recombination for equal mobilities : analytic solution without recombination . a) V = 
0.45, b) \q = 0.64, c) \ {) = 0.85, d) \ n =1.08 {for an ionization chamber of 1 cm electrode distance and a dose rate of 1 
R /sec , V = / gives a voltage of 4 . 35 volt ). 
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calculus which is straightforward giving in figure 4 only the results obtained for the space charge 
density at different applied voltages. Comparison with the curves of figure 3 shows that with 
recombination the maxima are very much lower and occur nearer to the electrode. Thus 
recombination displaces the depth of the maxima toward the center of the chamber and reduces 
their amplitudes. 

The effect is reflected in the behavior of the difference of space-charge densities u = p + - 
p_. For low values of the saturation ratio, u(x) increases nearly linearly with x. As i/i s increases, 
p+ increases in the cathode region and begins to be depleted in the anode region while the 
contrary happens to p_. This is shown in figure 5, again for various values of the saturation 
ratio, including recombination. 

Maxima of the space-charge density have been shown to occur for a = and a = 0.662. 
The situation is different for much higher values of a. It can be shown f.i. that for ex, — 1 
maxima do not occur any more. 



U 




FiGL'RK 5. Charge profiles- difference of space charge densities versus depth: computer solution for equal mobilities 
and a = .662 . A ) i/i s = 0.6,B)i/i s = 0. 75, C ji/i s = 0.99 (K = /*- = 1 .89XlO~ 4 m 2 l V sec). 
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5.2. Different Carrier Mobilities 

We have extended the computation of charge densities so as to include also the case of 
different carrier mobilities, figure 6 gives results obtained for the mobility values and 
recombination coefficient of Scott and Greening [5] and the saturation ratios of i/i s = 0.68, 0.77, 
and 0.98. Figure 6 shows that the curves for the less mobile positive carriers do not have 
maxima within the given interval of saturation ratios. The space charge density of the more 
mobile negative carriers does, however, exhibit a maximum for the lowest saturation ratio. 

For equal mobilities, symmetry demands that the net charge density becomes zero at x = 
1/2. Figure 7 shows that for different mobilities the zero is slightly displaced toward the anode, 
but still is practically independent of the saturation ratio. Mie has discussed this point [1], page 
884. His first approximation x Q = (k — 1)/(k 4- 1) is in good agreement with the computer results. 
Agreement becomes less satisfactory for his second approximation (eq 60), as can also be 
verified by inspection of figure 5, reference 4. 
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FlGURK 6. Charge profiles positive and negative charge densities versus depth. Computer solution for different 
mobilities (fji' + = 1.59X10 *m l \V sec, >jl'_ = 2.21 X 10 A m'W sec) and a = 0.7999. A) i/i s = 0.68, B) i/i s =0.77, C) i/i s = 
0.98. a) positive charge density; b) negative charge density. 
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Figikk 7. Charge profiles-net space charge density versus depth. Computer solution-parameters as in figure 6. 



6. Conclusion 

The present treatment shows that the positive and negative space-charge densities exhibit 
maxima which are due to the variation of the carrier velocities in the space-charge modulated 
electric field. 

Calculations were carried out on a Hewlett Packard 9100 A calculator, using the Runge- 
Kutta method of integration. 
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TABLE 1. List of quantities and symbols 



Quantity 



Symbol and 
dimension 



Reduced (nondimensional) 

quantity 



Electronic charge 

Dielectric constant 
Mobility of positive 

carriers 
Mobility of negative 

carriers 
Mobility ratio 
Recombination 

coefficient 
Rate of ionization 
Length 

Electrode distance 
Space charge density 
Difference of space 

charge densities 
( iurrent density 
Saturation current 

density 
Saturation ratio 
Electric field 

strength 
Voltage 
Applied voltage 



e = 1.602 x 10 n C 
e = 8.85 X 10 12 F/m 

/x+ (m fV sec) 

/Ut'_ (m \V sec) 



a (m /sec) 

q (C/m sec) 

x (m) 

d(m) 

()' (C/m 3 ) 



/" (A/rn ) 

/".(A/m 2 ) 
Hi. 

E (Vim) 

i" < n 
t'n(n 



AX_ = K 

« = a'eli'fjL + 

x = .v'/c/ 

p = p'(fxjeq) 112 

u = p + p 
i = i'lqd 

i.= i'Jqd = 1 

i = r/i' a 

E = dp J q) m E'ld 

v= (€ii J q) m rid 2 

V =(€fJLjq) ll2 VJ d 2 
/•> - E/i 
V*= V/i 2 



(Paper 79B1&2-423). 
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